EPJ manuscript No. 

(will be inserted by the editor) 



Classical Spin Liquid Properties of the Infinite-Component Spin 
Vector Model on a Fully Frustrated Two Dimensional Lattice 

Benjamin Canals^ and D. A. Garanin^ 

^ Laboratoire Louis Neel, CNRS, 25 avenue des Martyrs, Boite Postale 166, 38042 Grenoble Cedex 9, France 
^ Max-Planck-Institut fur Physik komplexer Systeme, Nothnitzer Strasse 38, D-01187 Dresden, Germany 

Received: date / Revised version: date 

Abstract. Thermodynamic quantities and correlation functions (CFs) of the classical antiferromagnet on 
the checkerboard lattice are studied for the exactly solvable infinite-component spin-vector model, D oo. 
In contrast to conventional two-dimensional magnets with continuous symmetry showing extended short- 
range order at distances smaller than the correlation length, r < oc exp(r*/T), correlations in the 
checkerboard-lattice model decay already at the scale of the lattice spacing due to the strong degeneracy 
of the ground state characterized by a macroscopic number of strongly fluctuating local degrees of freedom. 
At low temperatures, spin CFs decay as (SoSr) (x in the range ao ^ r ^ oc T~^^^, where ao is 
the lattice spacing. Analytical results for the principal thermodynamic quantities in our model are very 
similar with MC simulations, exact and analytical results for the classical Heisenberg model {D — 3) on 
the pyrochlore lattice. This shows that the ground state of the infinite-component spin vector model on 
the checkerboard lattice is a classical spin liquid. 
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1 introduction 

Frustrated magnets are particularly interesting as they of- 
ten involve unusual low temperature magnetic behaviors. 
Depending on the nature of interactions, the local con- 
nectivity, spin dimensionality, these systems stabilize un- 
conventional magnetic ground states such as non coUinear 
Neel orders, topological spin glasses, classical or quantum 
spin liquids, spin ices These observations are both ex- 
perimental and theoretical eventhough some of them, like 
topological spin glasses, are still conjectures. Among these 
ground states, two families have a striking property. Spin 
liquids, as well as spin ices, possess a residual entropy at 
zero temperature. This anomaly and the paradox it raises 
when considering the third principle of thermodynamics, 
are major points that people have tried to answer. 



spin correlations functions are exponentially decaying at 
finite temperatures, even sometimes at T = 0, in both 
classical and quantum cases. 



Recently, another lattice has received attention as it 
can be described as the two dimensionnal analog of the 
pyrochlore lattice : the checkerboard lattice. One very in- 
teresting point in this lattice is that its geometry mimics 
the local environment of each site of the pyrochlore lat- 
tice, but is much more simpler as it is two dimensional, and 
based on a scnaare lattice structure. This has allowed Lieb 
and Schupp ||] to establish exact results in the quantum 
case {S — 1/2) showing that the spectrum in this system 
is very peculiar as ground states are necesseraly singlets 
in finite size clusters. For > 1, it has been shown numer- 



Theoretical studies of frustrated antiferromagnets startedically M within an l/S expansion that the ground state 



half a century ago with the exact solution of the triangu- 
lar Ising antiferromagnet by Wannier . In the following 
years, many lattices have been identified where spin mod- 
els, from Ising symmetry {Z2) to Heisenberg symmetry 
(0(3) or SU(2)), have disordered ground states. Among 
them, two have attracted a lot of interest : the kagome 
lattice and the pyrochlore lattice The first one is a 
two dimensional arrangement of corner sharing triangles, 
while the latter is a three dimensional structure of corner 
sharing tetrahedra. When considering antiferromagnetic 
nearest neighbour interactions, they both display spin liq- 
uid like behaviors, well characterized theoretically. Spin- 



may have a local magnetization at zero temperature while 
for S = 1/2 it was proposed that the system orders in a 
valence bond solid ground state. A recent work Q derived 
a phase diagram of the ground states for varying S stating 
that this system is most likely ordered for all S. Exact di- 
agonalization of finite clusters |6| for S* = 1/2 reached the 
same conclusion and define the ground state as a 4-spin 
valence bond solid. A perturbative expansion also found 
a 4-spin valence bond solid [^| but with a different unit 
cell. There are now strong evidences that for small S this 
system should have long range order in a valence bond 
solid order parameter. 
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In this paper, we look for an exact solution on this 
lattice in the opposite limit of high spin dimensionality. 
Therefore, we study a classical Hamiltonian that corre- 
sponds to the generalisation of the Heisenberg model with 
I?-component spin vectors 



1 



and taking the limit D ^ oo. In this limit the problem 
becomes exactly solvable for all lattice dimensionalities, 

d, and the partition function of the system coincides 
with that of the spherical model. [ pl]Jl^ The D = oo 
model properly accounts for the profound role played, es- 
pecially in low dimensions, by the Goldstone or would 
be Goldstone modes. At the same time, the less signif- 
icant effects of the critical fluctuation coupling leading, 

e. g., to the quantitatively different nonclassical critical in- 
dices, die out in the limit D oo. Thus this model is 
a relatively simple yet a powerful tool for classical spin 
systems. It should not be mixed up with the iV-flavour 
generalization of the quantum S — 1/2 model |l3| i n the 
limit N — > oo, including its expansion. [[l4]Jl5| The 
A^-component nonlinear sigma model (see, e.g., Refs. 11161, 



as well as Ref. ||1 
extension of Eq. 



18[ for the 1 /N expansion) is a quantum 
in the long- wavelength region at low 
temperatures. Effective free energies for the n-component 
order parameter appear, instead of Eq. (^, in conven- 
tional theories of critical phenomena. Using them for the 
1/n expansion (see, e.g., Ref. is a matter of taste. 

While yielding the same results for the critical indices as 
the lattice-based 1/D expansion, [^0|j2l|p^ it misses the 
absolute values of the nonuniversal quantities. 

In the following, we give the solution of this model on 
the checkerboard lattice and show that its ground state is 
a classical spin liquid. Particularly, it is shown that many 
properties are similar to the pyrochlore lattice case de- 
spite the difference of lattice dimensionality. This strongly 
suggests that these frustrated systems are mainly driven 
by their local environment, i.e by their topological frus- 
tration, at least for D ^ oo. The rest of this article is 
organized as follows. In Sec. ||the structure of the checker- 
board lattice and its collective spin variables are described. 
In Sec. 1^ the formalism of the D = oo model is tailored 
for the checkerboard lattice. The diagrams of the classical 
spin diagram technique that do not disappear in the limit 
D oo are summed up. The general analytical expres- 
sions for the thermodynamic functions and spin CPs for all 
temperatures are obtained. In Sec. |^ the thermodynamic 
quantities of the checkerboard antiferromagnet are calcu- 
lated and compared with MC simulation results as well 
as exact and analytical results previously obtained on the 

r rochlore lattice in the whole temperature range. In Sec. 
the real space correlation functions are computed. We 
finally discuss our results in Sec. |^ and conclude. 

2 Lattice structure and the Hamiltonian 

Checkerboard lattice shown in Fig|| consists of corner- 
sharing tetrahedra. Each node of the corresponding Bra- 






Fig. 1. The checkerboard lattice. It can be pictured as a square 
lattice of tetraheda (top). Locally, it reproduces the same en- 
vironment as the three dimensional pyrochlore lattice. In our 
calculations, it has been described by a square lattice with a 
2-spins unit cell (bottom). 

vais lattice (i.e., each elementary tetrahedron in Pig.|l|) is 
numbered by i,j = 1,... ,N. Each site of the elemen- 
tary tetrahedron is labeled by the index / — 1,2. It is 
convenient to use the dimensionless units in which the in- 
teratomic distance equals 1/2 and hence the lattice period 
equals 1. The tetrahedra numbered hy i, j = I, . . . , N can 
be obtained from each other by the translations 

nuU + riyV, (2) 



J 



J 



where is the position of a site on the lattice, n„ and 
riy are integers, u and v are the elementary translation 



vectors (lattice periods), and 

u=(l,0), v = (0,l). 



(3) 



To facilitate the diagram summation in the next sec- 
tion, it is convenient to put the Hamiltonian (^ into a 
diagonal form. First, one goes to the Fourier representa- 
tion according to 



E 



s-e" 



(4) 



where the wave vector q belongs to the square Brillouin 
zone Bzkx,ky = [— tt, tt] x [— tt, tt] (see Pig. The Fourier- 
transformed Hamiltonian reads 



1 



2N 2^ ^q'S-q' 
H'q 



(5) 
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where the interaction matrix is given by 

with a — cos{qx), b = cos{qy), u — cos{qx/2) and v = 
cos(%/2). 

At the second stage, the Hamiltonian (|^) is finally di- 
agonalized to the form 

nq 

where — 2Jt'„(q) are the eigenvalues of the matrix 
taken with the negative sign, 

vi^l, j/2 = -(1 + cos((7j;) + cos((7y)). (8) 

The diagonalizing transformation has the explicit form 

C/-/(q)<f/,,„.(q) = Vq",S„„S (9) 

where the summation over the repeated indices is implied 
and tj is the real unitary matrix, tj~^ — U"^ , i.e., J7~/ = 
Uin- The columns of the matrix U are the two normalized 
eigenvectors J7„ = (?7i„, U2n) of the interaction matrix V: 

/ l+b fl+a \ , , 

V 2 + a + \ 2uv } 

The normalized eigenvectors satisfy the requirements of 
orthogonality and completeness, respectively: 

C^in(q)t^i«'(q) = ^rm', C^/n (q) C/;'„ (q) = . (11) 

The Fourier components of the spins Sq and the collective 
spin variables are related by 

The largest dispersionless eigenvalue v\ of the interac- 
tion matrix [see Eq. (H)] manifests frustration in the sys- 
tem which precludes an extended short-range order even 
in the limit T — > 0. Independence of v\ of q signals that 
1/2 of all spin degrees of freedom are local and can rotate 
freely. The other eigenvalue satisfy 

i.2(q)^-3 + <zV2 (13) 

at small wave vectors, <f = + Qy <C 1. Within our 
description of the checkerboard lattice, the would be usual 
Goldstone mode corresponds to 1^2- At the corners of the 
Brillouin zone (±7r, ±7r), this mode becomes degenerate 
with the flat one and defines low energy excitations. 

l^2{ci) = I - qV2 (14) 




Fig. 2. Reduced eigenvalues of the interaction matrix, i^n(q) = 
V'q /(2J) of Eq. (Ql, plotted over the Brillouin zone. The flat 
band signals that half of the degrees of freedom in this system 
can rotate freely. 

where q — (±7r, ±7r) -|- q, \q\ << 1. It is the equiva- 
lent of the usual Goldstone mode that is present in low- 
dimensional magnets with a continuous symmetry. Thus, 
it is this mode that will be associated with a correlation 
length, as we will see in section^. It is remarkable that this 
mode contains all information about dispersive excitations 
present in this model. If we had choosen a four sites unit 
cell [|4) instead of a two sites one, we would have obtained 
a similar description as the one done on the pyrochlore an- 
tiferromagnet [^. Then, an "optical" mode would have 
been recovered, i.e a mode with a flnite energy gap for all 
wave vector; the usual Goldstone mode would have been 
degenerate with the flat one at g = and the lowest not 
dispersive branch, would have been two times degenerate. 
This clearly shows the similarity of these two lattices as it 
is expected because of the exact equivalence of the local 
environment of each site. The q dependences of the eigen- 
values Vn over the whole Brillouin zone are shown in Fig. 

i 

In the next section the equations describing spin correla- 
tion functions of the classical checkerboard antifcrromag- 
net in the large- 1? limit will be obtained with the help of 
the classical spin diagram technique. The readers who are 
not interested in details can skip to Eq. ( [l^ or directly to 
Sec. |. 

3 Classical spin diagram technique and the 
large-D limit 

The exact equations for spin correlation functions in the 
limit D ^ 00, as well as the 1/D corrections, can be the 
most conveniently obtained with the help of the classical 
spin diagram technique. [ p6|j2^ , ^ A complete description 
of the technique applied to a non bipartite lattice can be 
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found in Ref. ||2^. We only give here an outline of the 
calculations. 

Our goal is to compute spin-spin CF's 



s"'(q) 



which are related to CF's of the a variables 



„n / \ ^ i n n 



-q) = ]^(^q^-c 



(15) 



(16) 



(17) 



following from Eq. (^. 

The analytical expression for the a CF in the SCGA 
has the Ornstein-Zernike form 



through the relation 



s"'(q) = t/,„(q)C/r„(q)a"(q) 



a"(q) 



DA 



1 - Am 



In the large-D limit the expression of A simplifies 
is given by 



A = 



1 



D 1 + yr+siTD 

Here the dispersion L is given by the formula 

dq (/?Vq")' 



L 



A 



2-2! 



E/" dc 
. "7(^ 



(^■kY 1 - AfiV^ 



(18) 



and 



(19) 



(20) 



The summation (1/A^) X]q ■ • ■ is replaced by the integra- 
tion over the Brillouin zone, vq is the unit cell volume, 
and d is the spatial dimensionality. For the square lattice 
we have = 1. The expression for L can be simplified to 



L = 



P~ 1 
'2A~'' 



(21) 



where Pn is the lattice Green function associated with the 
eigenvalue n: 



Pn = vo 



dq 



1 



i2TTY 1 - APV^ 



(22) 



Now one can eliminate L from Eqs. (|T^) and (|2l|), which 
yields the basic equation of the large-D model, 



DAP ^ 1. 



(23) 



sites r. Indeed, calculating the spin autocorrelation func- 
tion in the form symmetrized over sublattices with the 
help of Eqs. (|l^, ([u]), and (^, one obtains 



'jj^lE'-M-'^'ifi^ (24) 



Wo 



That is, the spin-normalization condition is automatically 
satisfied in our theory by virtue of Eq. (^3|). After A has 
been found from this equation, the spin CFs are readily 
given by Eqs. (^ and (y). 

To avoid possible confusions, we should mention that 
in the paper of Reimers, Ref. where Eq. (^8|) with 
the bare cumulant A= 1/D has been obtained, the theo- 
retical approach has been called "Gaussian approximation 
(GA)". This term taken from the conventional theory of 
phase transitions based on the Landau free-energy func- 
tional implies that the Gaussian fluctuations of the order 
parameter are considered. In the microscopic language, 
this merely means calculating correlation functions of fluc- 
tuating spins after applying the MFA. Such an approach 
is known to be inconsistent, since correlations are taken 
into account after they had been neglected. As a result, 
for the checkerboard lattice one obtains a phase transition 
at the temperature = T^^^ = 2J/D but immediately 
finds that the approach breaks down below Tc because of 
the infinitely strong fluctuations. In contrast to this MFA- 
based approach, the self-consistent Gaussian approxima- 
tion used here allows, additionally, to the Gaussian fluc- 
tuations of the molecular field, which renormalize A and 
lead to the absence of a phase transition for this class of 
systems. The SCGA is, in a sense, a "double-Gaussian" 
approximation; one concerning fluctuations of the order 
parameter and the other one describing fluctuations of 
the molecular field. 

To close this section, let us work out the expressions 
for the energy and the susceptibility of the checkerboard 
antiferromagnet. For the energy of the whole system, using 
Eqs. (0) and (|l|), as well as the equivalence of all spin 
components, one obtains 



(2^)-^ 1 



(q). (25) 



To obtain the energy pro spin U, one should divide this 
expression by 2A^. With the use of Eq. (18), the latter can 



be expressed through the lattice Green's function P of Eq. 
. . (eTI) ; then with the help of Eq. (E3) it can be put into the 

This nonlinear equation determining yl as a function of^ jinal form 

temperature differs from those considered earlier 2^p9| , p7[ f 
by a more complicated form of P reflecting the lattice 



structure. The form of this equation is similar to that ap- 
pearing in the theory of the usual spherical model. ||ll| , ^ 
The meanings of both equations are, however, different. 
Whereas in the standard spherical model a similar equa- 
tion account for the pretty unphysical global spin con- 
straint, Eq. ( p3|) here is, in fact, the normalization condi- 
tion (s^) = 1 for the spin vectors on each of the lattice 



D- ^ 
A 



(26) 



The susceptibility pro spin symmetrized over sublattices 
can be expressed through the spin CFs as 



1 



Xq = 



2DT 



(27) 



Benjamin Canals, D. A. Garanin: Title Suppressed Due to Excessive Length 



5 



With the use of Eq. ( pT| ) this can be rewritten in the for: 
n ; 

(21 

where the diagonalized CPs are given by Eq. (|l^). Fro: 
Eq. (|o|) it follows that in the hmit q — > one has Wi — 
and W2 — \/2. Thus the homogeneous susceptibility x 
Xo simplifies to 

X=^'t'{0). (2! 

As we shall see in the next section, disappearance of tl 
terms with n = 1 from this formula ensures the nondive 
gence of the homogeneous susceptibility of the checke 
board antiferromagnet in the limit T — > 0. The situatic 
for q ^ is much more intricate and it will be considered 
below in relation to the neutron scattering cross section. 




Fig. 3. Temperature dependence of the gap parameter G for 
the checkerboard antiferromagnet. 



4 Thermodynamics of the checkerboard 
antiferromagnet 

To put the results obtained above into the form explicitly 
well behaved in the large-Z) limit and allowing a direct 
comparison with the results obtained by other methods 
for systems with finite values of Z?, it is convenient to use 
the mean-field transition temperature T^^^ = 2J/D as 
the energy scale. With this choice, one can introduce the 
reduced temperature 9 and the so-called gap parameter G 
according to 



T 



~ yMFA' 

In these terms, Eq. ( p3| ) rewrites as 
OGPiG) = 1 



G.^A. 



(30) 



(31) 



and determines G as function of 9. Here P{G) is defined 
by Eq. (^l|), where 



Pi 



1 



1-G" 



P2 vo 



1 



(27r)'il-Gi.2(q)' 



(32) 



The (T CPs of Eq. (|l8|), which are proportional to the in- 
tegrands of P„, can be rewritten in the form 



a"(q) = 



9G 



1-Gi.„(q)' 



(33) 



Further, it is convenient to consider the reduced energy 
pro spin defined by 



U = U/\Uo\ 



Uo 



-J, 



(34) 



where Uq is the energy pro spin at zero temperature. With 
the help of Eq. (Bq) U can be written as 



U 



1/G. 



(35) 



The homogeneous susceptibility x of Eq. ( |29D can be rewrit- 
ten with the help of Eq. (p^ in the reduced form 



G 



1 + 3G 



(36) 



The sense of calling G the "gap parameter" is clear 
from Eq. (33). If G = 1, then the gap in correlation 
functions closes: turns to infinity, and diverges at 
q — !■ (±7r, ±7r). For nonordering models, it happens only 
in the limit 6^0, however. The solution of Eq. ( ^ ) sat- 
isfies G < 1 and goes to zero at high temperatures. If 
<C 1, the function P is dominated by Pi = 1/(1 — G). 
The ensuing asymptotic form of the gap parameter at low 
temperatures reads 



G^ 1 



< 1. 



(37) 



At high temperatures, Eq. ( |3l|) requires small values of 
G. Here, the limiting form of P can be shown to be P = 
1 -I- (3/2)G^. The corresponding asymptote of G has the 
form 



G = 



1 



3 
2^ 



> 1. 



(38) 



The numerically calculated temperature dependence of G 
is shown in Fig. |. Note that in the MFA one has G = 1 /9 
which attains the value 1 at 6* = 1. 

The temperature dependence of the reduced energy of 
Eq. (|3^) is shown in Fig. ^. Its asymptotic forms following 
from Eqs. (p^) and (pq) are given by 



U' 



-3/{29), 
-1 + 9/2, 



> 1 
< 1. 



(39) 



This implies the reduced heat capacity G — dU jdB equal 
to 1/2 at low temperatures. To compare with the MC 
simulation data of Ref. pit] for the heat capacity of the 
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-0.8 



Fig. 4. Temperature dependence of the reduced energy for th( 
checkerboard antiferromagnet. 



0.80 



C = dU/dT 




e = T/T, 



Fig. 5. Temperature dependence of the heat capacity of the 
checkerboard antiferromagnet. The MC results of Ref. for 
the Heisenberg model {D — 3) on the pyrochlore lattice are 
represented by circles. Black line corresponds to the exact so- 
lution of the infinite component heisenberg antiferromagnet on 
the pyrochlore lattice . Dashed line corresponds to the exact 
solution of the infinite component heisenberg antiferromagnet 
on the checkerboard lattice. 



Heisenberg model we will use, instead of C, the true heat 
capacity C = dU/dT = [D/2)C [see Eqs. @ and (||)], 
which in our approach tends to D/4 3/4 at low temper- 
atures. The temperature dependance of the heat capacity, 
compared with previous MC results of the Heisenberg an- 
tiferromagnet on the pyrochlore lattice |^ as well as exact 
solution of the infinite component classical antiferromag- 
net on the pyrochlore lattice |Q are shown in Fig. p|. We 
see that the behavior of these two different models pl| , p3[ 
on the pyrochlore lattice are very similar to our results on 
the checkerboard lattice, thus confirming the analogy of 
these two structures despite their different space dimen- 
sionality. 



0.30 



0.20 



0.10 



0.00, 




e = t/t: 



Fig. 6. Temperature dependence of the reduced uniform sus- 
ceptibility of the checkerboard antiferromagnet. Asymptotes 
for high and low temperatures are given, for both the infinite 
component model and the mean field case. 



Using Eq. ^ we compute the low and high tempera- 
ture asymptotic behavior of the reduced susceptibility x 



X 



f/4- 61/32, 
(l/^?)(l-3/0). 



< f 
> 1. 



(40) 



Its dependence at all temperatures is reported in Fig. g. 
Here also, we recover a behavior very similar to the one 
of the pyrochlore lattice where MC simulations and ana- 
lytical results have been obtained for the Heisenberg an- 
tiferromagnet as well as exact results for the infinite 
component classical antiferromagnet p3[ |. The similarity 
between our results and the local approximation [33| de- 
velopped in Ref. suggest that the checkerboard an- 
tiferromagnet has its thermodynamics governed by local 
correlations, as it is expected in a classical spin liquid. 



5 Real-space correlation functions 

The low-temperature behavior of the a correlation func- 
tions of Eq. (^3|) is dominated by their asymptotic form at 
small wave vectors, where the would be Goldstone mode 
is defined, i.e at q =_(±7r, ±7r)_f g, |g| << 1. According to 
Eqs. (|), we find, and (||), by 



2k2 



(41) 



where the quantity = 6 in defines the correlation 
length 



(42) 
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Appearance of this length parameter implies that the real- 
space spin CFs defined, according to Eqs. (|l7|) and (0), by 



Ji' 



(2^ 



C/;„(q)(7/'n(q)a"(q) (43) 



decay exponentially at large distances at nonzero temper- 
atures, fn contrast to conventional lattices, divergence of 
at ^ does not lead here to an extended short- 
range order, i.e., to strong correlation at distances r < £^c- 
The zero-temperature CFs are purely geometrical quanti- 
ties which are dominated by and have the form 



2vQ 



dq 



.e'i<-\~-\)Un{ci)Ui.M) (44) 



It is convenient to enumerate CFs by the numbers n„ 
and defined by Eq. (||). Thus s"^ „^ is the correlation 
function of the I sublattice spin of the "central" unit-cell 
(0, 0) with the sublattice spin of the unit-cell translated 
by (n„,n«). Let us calculate the CFs with I = V = 1 &i 
large distances along the horizontal line in Fig. [l], at small 
but non zero temperature. We use the asymptotic form of 
the matrix [/i,i(q) for q ~ (±7r,±7r) -|- q, \q \ « 1, as well 
as the asymptotic form of CFs in Eq. ^ This allows to 
write this correlation function as 



2vQ 



dq 

(2^ 



-1)V"«==(^ + ||V 



g2. 



(45) 



Simplifying this expression and taking into account that 
the integral of the cosine over the whole Brillouin zone is 
zero, one arrives at the asymptotic form 



(-1)^ 



(-1)"k2 



/27r (Kn)3/2 



, Kn » 1, 



(46) 



where K^(x) is the Macdonald (modified Bessel) function. 
The strong decrease of this correlation function with dis- 
tance even at T = (k = 0) is not surprising, since our 
solution spans the whole highly degenerate ground-state 
manifold and this degeneracy is not lifted in the limit 
D — > oo. At small but non zero temperature, spin-spin 
correlations are exponentially decaying as if interactions 
were renormalized to zero and drive this system to a para- 
magnetic fix point for finite temperatures. Let us calcu- 
late more generally the CFs with I = 1,2 and Z' = 1, 2 at 
large distances along the horizontal line in Fig. |], at zero 
temperature. Following the same previous procedure, we 
obtain 



oil 



„22 



7r(n + 1/2)2 



(47) 



where n is in unit of the inter-cell distance. We note that 
despite the divergence of the correlation length when 



n+1,0 




0.01. 



0.001 r 



0.0001 4 



n+1 (inter-cell spacing) 

Fig. 7. Real-space correlation functions sJ,^o at T = calcu- 
lated from Eq. along an horizontal line. The distance unit 
is the interatomic spacing. The asymptote given by Eq. (^6|) is 
shown by the dashed lines. 



lowering the temperature, the model do not order, which 
clearly defines this system as a classical two dimensionnal 
spin liquid. 



6 Discussion 

In the main part of this article, we have presented in de- 
tail the exact solution for the D = oo component classi- 
cal antiferromagnet on the checkerboard lattice. The so- 
lution does not show ordering at any temperature due to 
the strong degeneracy of the ground state, and thermo- 
dynamic functions behave smoothly. In contrast to con- 
ventional two-dimensional magnets, there is no extended 
short-range order at low temperature, and T = is not 
a critical point of the system. Although the correlation 
associated to the would be Goldstone mode diverges as 
C!C T"^/^, the power law decay (soSr) « 1/r^ at zero 
temperature of the spin correlation functions leads to the 
loss of correlations at the scale of the interatomic dis- 
tance. All these properties characterize this model as a 
two-dimensional classical spin liquid. As these results are 
obtained within the D — oo component classical antifer- 
romagnet, it is expected that it mimics the large S value 
Heisenberg antiferromagnet on the checkerboard lattice. 
Surprisingly, there are results |^,^ suggesting that for 
large S, the Heisenberg antiferromagnet should be ordered 
on this lattice. At this stage, the discrepancy between the 
different approaches is not clear. A possible explanation , 
as reported in Ref. "is that for large S, there should be 
quantum order by disorder (as supported in ||^) but this 
order by disorder comes with an energy scale of J/S rather 
than J. Therefore, in the classical limit, the temperature 
at which the correlation length grows exponentially goes 
to zero. In other words, the limits T and 5' ^ oo (or 
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D oo) should not commute." Nevertheless, it should be 
stressed that even if there is a discrepancy at T 0, for 
finite T, we expect the D = oo case to correctly describe 
the checkerboard Heisenberg antiferromagnet as any crit- 
ical behavior is expected to appear only in the zero tem- 
perature neighbourhood. A similar problem has been en- 
countered within the present formalism applied to kagome 
antiferromagnet p8| . Here again, finite temperature prop- 
erties of the infinite component spin vector model have 
been found to be close to the D = 3 Heisenberg case. How- 
ever at T — > 0, the method misses the "order by disorder" 
phenomenon which selects the coplanar spin manifold. As 
a consequence, the _D = 3 zero temperature specific heat is 
much closer to 1, C = ll/12fcB, than it would be if really 
one third of the degrees of freedom were still fiuctuating 
(as suggested by a mean field analysis). Wether 1/D cor- 
rections could resolve this discrepancy is still a work to be 
done. 

Despite the dimensionality of the checkerboard lattice 
{d = 2), there are many similarities with previously ob- 
tained results on the pyrochlore lattice (c? = 3). The spe- 
cific heat has in particular, the same value at zero tem- 
perature [^3[|3l|] for D ~ 3. Susceptibility is also very 
similar with the one of the Heisenberg antiferromagnet as 
well as the infinite component classical antiferromagnet 
on the pyrochlore lattice. MC results, as well as exact and 
analytical results are very close to the one obtained here 
1^,^. This is not surprising as the spin liquid behav- 
ior implies that the lattice is mainly described by local 
fluctuations. In fact, it is true at least for finite tempera- 
ture thermodynamics. Why the present approach can lead 
to correct results for the pyrochlore case and not for the 
checkerboard case at very low temperatures has not been 
addressed in the present work. We can only note that the 
common point between the two lattices is the local con- 
nectivity which is clearly well taken into account within 
the D ^ oo formalism. Wether dimensionality is relevant 
is not clear. At this stage, it can be noted that for the 
kagome case, the D ^ oo formalism has missed the re- 
duced specific heat mechanism (see Ref. ||2^,Q) while for 
the pyrochlore case, it has been reproduced (see Ref. psj) 
Therefore, discrepencies with D = 3 descriptions could be 
ascribed to dimensionality although other subtle phenom- 
ena that are not taken into account when D oo cannot 
be ruled out without studying the same model when in- 
cluding 1/D corrections. 

To conclude, this classical approach indicates that for 
infinite component spin value, the checkerboard antiferro- 
magnet should be disordered and behave as a spin liquid. 
We give some arguments why the present approach is at 
variance with results of Ref. It is interesting to note 
that for small spin value, the checkerboard antiferromag- 
net is expected to order ||,^,^ . Therefore the checkerboard 
antiferromagnet cannot be compared to the pyrochlore an- 
tiferromagnet even if its geometry defines it as its two di- 
mensional analog. Even if they have the same unit cell, 
it is probably the global geometry of the underlying Bra- 
vais lattice that drives the physics and not only the local 
connectivity. 
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